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ABSTRACT 

We analyse simultaneous UBVR quiescent light curves of the cataclysmic variable 
V2051 Oph using the Physical Parameter Eclipse Mapping method in order to map 
the gas temperature and surface density of the disc for the first time. The disc appears 
optically thick in the central regions and gradually becomes optically thin towards the 
disc edge or shows a more and more dominating temperature inversion in the disc 
chromosphere. The gas temperatures in the disc range from about 13 500 K near the 
white dwarf to about 6 000 K at the disc edge. The intermediate part of the disc has 
temperatures of 9 000 K to 6 500 K. 

The quiescent disc (chromosphere) shows a prominent bright region with temper- 
atures of 10 500 K around the impact region of the stream from the secondary with an 
extension towards smaller azimuths. The disc has a size of 0.53 ± Q.03R L1 and a mass 
accretion rate of between M. = 10 15 gs _1 to 10 17 gs _1 . The light curves must include 
an uneclipsed component, a hot chromosphere and/or a disc wind. 

The PPEM method allows us to determine a new distance of 146 ± 20 pc, com- 
patible with previous rough estimates. For the white dwarf we then reconstruct a 
temperature of 19 600 K, if the lower hemisphere of the white dwarf is occulted by the 
disc. 

We suggest that the accretion disc is a sandwich of a cool, optically thick central 
disc with hot chromospheric layers on both sides as was suggested for HT Cas. This 
chromosphere is the origin of the emission lines. 

We find that although V2051 Oph is very similar to the SU UMa type dwarf novae 
HT Cas, OY Car and Z Cha, there must be a substantial difference in order to explain 
its unique light curve. The reason for the difference could either be a higher mass 
transfer rate caused by the more massive secondary and/or a small but significant 
magnetic field of the white dwarf, just strong enough to dissrupt the innermost disc. 

Key words: binaries: eclipsing - novae, cataclysmic variables - accretion, accretion 
discs - stars: V2051 Oph 



1 INTRODUCTION 

The eclipsing cataclysmic variable (CV) V2051 Oph was dis- 
covered by Sanduleak (1972). Since then, the classification 
as a CV was quite certain, but not, of which subtype. Bond 
& Wagner (1977) suggested it might be of AM Her type and 
this idea was supported by Warner & O'Donoghue (1987), 
hereafter WO, refining the classification to a low-field polar. 
Other studies are based on the assumption that V2051 Oph 
is a dwarf nova (Watts et al. 1986, Baptista et al. 1998) 

* Send offprint requests to: S. Vrielmann 



with infrequent outbursts and low states (Warner & Crop- 
per 1983, Baptista et al. 1998). V2051 Oph shows typical 
features of a high inclination dwarf nova exhibiting an ac- 
cretion disc, like double peaked emission lines. 

The mystery of V2051 Oph seemed to be solved with 
the detection of a super outburst during which superhumps 
were detected (Kiyota & Kato 1998). Since then, the sys- 
tem was classified as an SU UMa type dwarf nova. There- 
fore, we assume this system to exhibit an accretion disc. 
This is further supported by Baptista et al. (1998) who do 
not find typical features of a polar in their light curves of 
V2051 Oph, like wanderings of ingress and egress phases or 
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Table 1. The quiescent light curves of V2051 Oph. 



Date 


UT(start) 


Eclipse No. 


14/6/83 


6:51:06 


1,2 


15/6/83 


6:35:29 


17 


16/6/83 


5:38:05 


32,33 


17/6/83 


5:28:07 


48 


18/6/83 


5:02:52 


64,66 



Table 2. The eclipse timings of our data. The cycle number and 
observed minus calculated times are given with respect to Bap- 
tista et al.'s ephemeris. 



large amplitude pulsations. However, in order to explain the 
complex behaviour of V2051 Oph, Warner (1996) suggests 
that this system is a polaroid, an intermediate polar with a 
synchronized primary. 

With 1.5 it has a rather short orbital period, very sim- 
ilar to HT Cas, OY Car and Z Cha and in fact it is often 
compared to these. However, the accretion disc seems to be 
highly variable, as apparent from the variable eclipse pro- 
files. Sometimes the disc seems to be only small, hot and 
sharply defined, while at other times the disc appears to be 
more extended (Warner & Cropper 1983). Furthermore, at 
times it shows a hump in the light curve with variable tim- 
ing (before or after eclipse) and variable occurance (it can 
occur at one eclipse and disappear at the next) (WO). At 
other times, no such hump is observed (e.g. this paper). 

As V2051 Oph is an eclipsing system, the system pa- 
rameters are quite well known (e.g. Baptista et al. 1998). 
However, the lack of traces of the secondary even in the 
infra-red prevents pinning down the distance. The estimate 
of Watts et al. (1986) of a distance between 90 to 150 pc is 
based on the assumption that the quiescent disc emits in an 
optically thick fashion. Berriman, Kenyon & Bailey (1986) 
derive values between 120 and 170 pc, allowing for optically 
thin material in the disc. 

One of our aims was to pin down the distance with the 
help of the Physical Parameter Eclipse Mapping (PPEM, 
Vrielmann, Home & Hessman 1999, VHH) method as ex- 
plained Vrielmann (2001) and briefly described in Sec- 
tion 



3.3. The PPEM is a Maximum Entropy method (MEM, 
Skilling & Bryan 1984) designed to analyse accretion discs 
in CVs: by analysing the shape of the eclipse light curve in 
various colours, we can retrieve information about the phys- 
ical properties within the accretion disc. For this, we assume 
a certain spectral model for the disc emission. The applica- 
tion of this method to quiescent data is shown in Section |5| 
The last Section H summarizes and discusses the results. 



2 THE DATA 

The data set was taken by R. Stiening between 14 and 18 
June, 1983 with the Mt. Lemon/USA 1.5m telescope. At- 
tached to the telescope was the SLAC (Stiening) photome- 
ter allowing the detection of light in the four filters UBVR 
simultaneously. The data set was calibrated and reduced by 
Keith Home and consists of eight eclipses with a time resolu- 
tion of about 1 sec. A log of the data is shown in Table [lj The 
procedure of the observations is similar to that of HT Cas 
by Home, Wood & Stiening (1991). 

For the analysis of the quiescent data we assume a flat 
disc geometry. This involves adjusting the out-of-eclipse lev- 



Cycle 


HJD 


(O-C) 




(2445000.+) 


(cycles) 


36103 


499.81007 


-0.0123 


36104 


499.87249 


-0.0125 


36119 


500.80902 


-0.0107 


36135 


501.80779 


-0.0119 


36150 


502.74424 


-0.0114 


36166 


503.74321 


-0.0094 


36168 


503.86801 


-0.0103 



els to the same level before and after eclipse by fitting a low 
order polynomial to the out-of-eclipse data. 



3 THE PPEM ANALYSIS 

A preliminary analysis of the data set was shown in Vriel- 
mann (1999). We arrive at some different results in the 
present study because of redetermination of the ephemeris 
(using the 1983 data alone and exclusion of two faulty R- 
band eclipses) and the correction of an unfortunate error in 
the computer programme which changed the results signif- 
icantly. A slight difference to the result presented in Vriel- 
mann (2001) is due to a careful reexamination of the original 
data. This, however, did not change the results significantly. 



3.1 The light curves 

Fig. ^ and ^| show the individual UBVR light curves of 
V2051 Oph at the full time resolution of Is, showing the 
eclipse of the white dwarf and accretion disc by the sec- 
ondary. The gaps are due to calibration measurements, un- 
fortunately often during mid-eclipse. However, since the 
ingress and egress contain most information on the disc emis- 
sivity, this choice for calibration measurements is a compro- 
mise one can live with. The R-band light curves no. 1 and 
no. 2 were excluded from the analysis due to missing or 
faulty sky measurements. 

The phases were calculated according to a newly deter- 
mined ephemeris of 

HJD(eclipse) = 2445499.81007 + 0.062430039 x E (1) 
BJD(eclipse) = 2445499.81002 + 0.062430039 x E (2) 

using 26 out of the 30 light curves displayed in Fig. [j] 
and ^ (eclipse no. 32 was excluded, because of the incom- 
plete coverage of the eclipse ingress) . It was derived from the 
mean of the steepest ingress and egress variation (at phases 
±0.033) which can be attributed to the eclipse of the white 
dwarf by the secondary. Table |^ gives the (O-C) values with 
respect to Baptista et al.'s (1998) ephemeris. Our (O-C) val- 
ues are somewhat large (although not unlike the scatter in 
other measurements) and not in agreement with Baptista's 
(2001, private communication) cyclic ephemeris. However, 
special care has been taken in deriving the mid-eclipse times 
(e.g. including leap seconds) and the error in our mid-eclipse 
times is certainly much smaller than the (O-C) values given. 

While the white dwarf eclipse is relatively stable (see 
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Figure 1. The individual U-band and B-band light curves. The numbers on the right give the eclipse numbers according 
The eclipse light curves were offset by 5 mjy. 



to Table 



Table 3. The white dwarf ingress and egress times and phases. 



HJD ingress 


ingress phase 


HJD egress 


egress phase 


(2445000.+) 




(2445000.+) 




499.80800 


-0.033189 


499.81214 


0.033125 


499.87044 


-0.033030 


499.87454 


0.032644 


500.80696 


-0.031918 


500.81108 


0.034076 


501.80574 


-0.033530 


501.80984 


0.032144 


502.74221 


-0.033219 


502.74628 


0.031974 


503.74106 


-0.033709 


503.74536 


0.035168 


503.86594 


-0.033390 


503.87007 


0.032764 




-0.03314 




0.03313 



Tab. ^|), the remainder of the light curve is quite variable, 
showing flickering and flares which cause the accretion disc 
ingress and egress profile to vary quite significantly. Dur- 
ing white dwarf eclipse the flares are less strong, though 
not completely absent, indicating that the region where the 
flares originate is not only confined to the immediate vicin- 
ity of the white dwarf, but also extends somewhat into the 
inner disc. 

The variation in the out-of-eclipse light curve is similar 
in all filters, e.g. flares occur in light curves no. 32 (phase 0.09 
and 0.28) or no. 66 (phase 0.12). In the V-band and R-band 



light curves the flares appear slightly stronger, however, not 
as extreme as in HT Cas (Home et al. 1991), where a strong 
flare present in the R-band appears nearly absent in U-band. 

It is interesting to note that even though we see strong 
variability of the eclipse light curves in form of flickering and 
flares, we do not see as strong and broad humps as in some of 
WO's light curves. Even though our data show a lower S/N 
ratio such strong humps cannot be hidden, but their absence 
reflect a seasonally changing behaviour of V2051 Oph. Since 
WO present their data in counts/sec without comparison 
with a star of constant flux, it is difficult to compare their 
data with ours. 

WO also claim not to see the steep ingress and egress 
features of the white dwarf, while we can certainly distin- 
guish quite well the white dwarf eclipse from the accretion 
disc eclipse in most individual light curves. However, a close 
look at their eclipse light curves shows in many cases a char- 
acteristic kink in the lower part of the steep ingress profile, 
possibly the second contact phase of the white dwarf, and 
in many cases a corresponding kink (third contact phase of 
the white dwarf) in the egress profile. 

In spite of the varying eclipse profile, we averaged the 
light curves in order to reduce the flickering which would 
otherwise introduce artefacts into the PPEM reconstruc- 
tions. To enforce the fit of the eclipse ingress, egress and 
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the eclipse depth we chose a phase resolution twice as high 
within the phase range —0.1 to 0.1 as outside. The total 
phase range from —0.15 to 0.15 excludes some flares, how- 
ever not all, noticable especially at phases —0.135 and 0.13. 
The large variation of the out-of-eclipse light curve levels is 
reflected in the error bars, which are therefore smaller during 
eclipse than outside the eclipse profile. 

The averaged and phase-binned light curves are dis- 
played in Fig. pi We fitted a low-order polynominal to the 
out-of-eclipse light curves and divided the full light curves 
by the polinomia in order to reduce any variation caused by 
an orbital hump, though it was not very obvious in the av- 
eraged light curve. An orbital hump, as often seen in dwarf 
nova light curves, is usually caused by anisotropic emission 
from the disc edge, in most cases maximally seen just pro- 
ceeding the eclipse. Since we used a geometrically thin ac- 
cretion disc geometry, we would not be able to account for 
such an effect. 

In comparison with the white dwarf ingress and egress 
profile, the disc eclipse is very asymmetric. Within the white 
dwarf eclipse (between phases —0.033 and 0.033) the disc 
ingress is rather V shaped, indicating a flat intensity pro- 
file, while the disc egress is U shaped, indicative of a steep 
intensity profile. This is apparent in all filters, though less 
pronounced in the noisier R-band light curve. Furthermore, 



the flux at white dwarf ingress is higher than at white dwarf 
egress. All these features are indicative of a brighter leading 
lune of the disc, i.e. a bright spot where the ballistic stream 
from the secondary hits the accretion disc. 



3.2 The PPEM method and the model 

The PPEM method is based on the eclipse mapping (EM) 
technique first introduced by Home (1985). In EM one takes 
advantage of the spatial information of the intensity distri- 
bution of the accretion disc hidden in the eclipse profile. By 
fitting the eclipse light curve of high inclination systems, 
caused by the occultation of the accretion disc by the sec- 
ondary star, we reconstruct the intensity distribution using a 
maximum entropy method (MEM) . The MEM algorithm al- 
lows one to choose the simplest solution still compatible with 
the data in this otherwise ill-conditioned back projection 
problem. Further details of the EM method can be found in 
Home (1985) or Baptista & Steiner (1993). 

As described in VHH, the further developed PPEM 
method is designed to reconstruct parameter distributions, 
like temperature T and surface density E, rather than inten- 
sity distributions of the accretion disc. Hereby, the assump- 
tion of a spectral model is necessary, relating the parameters 
to be mapped (e.g. T, E) to the radiated intensity in a given 
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UBVR light curves of V2051 Oph, x = 2.5 



Entropy vs. Distance 




Orbital Phase 

Figure 3. The averaged quiescent light curves in UBVR together 
with the fits. The ingress and egress phases of the white dwarf 
(determined as steepest ingress and egress at phases —0.033 and 
0.033) are drawn as vertical dashed lines. The data were offset by 
6 mjy each with dotted horizontal lines indicating the zero-level 
for each light curve and solid lines indicating the phase constant 
uncclipscd flux (for the V-band light curve they coincide) . Resid- 
uals are displayed as solid lines in relation to the zero-level lines. 



filter I v = f(T, E). One has to be cautious, however, as the 
parameters are not necessarily everywhere well defined: only 
after the reliability of the reconstructed parameters is estab- 
lished can the resulting maps be interpreted. 

We used the two-parameter spectral model described in 
VHH, to map the temperature and surface density of the 
disc material. Though it is a very simple, pure hydrogen 
model including only bound-free and free-free H and H~ 
emission, it is still very useful, because it allows us to distin- 
guish between regions of the disc that emit like a black body 
and those which deviate from it. In the simplest interpreta- 
tion, the latter are optically thin regions and presumably 
the site of line emission (Watts et al. 1986, Baptista et al. 
1998). However, real disc might show a Balmer Jump in ab- 
sorption due to decreasing temperature with disc height or 
a temperature inversion in the disc photosphere (i.e. a disc 
chromosphere with line emission). If such a more compli- 
cated scenario is present in the disc, our simple solution will 
fail. 

For the white dwarf emission we used white dwarf spec- 
tra. As described in VHH, we assumed the white dwarf to 
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Figure 4. The entropy of the reconstructions (as a measure of 
smoothness) as a function of the assumed distance. The dashed 
line is a parabolic fit to the data. It peaks at a distance of 146 pc. 



be spherical in terms of occultation by the secondary and 
occultation of the accretion disc by the white dwarf. Oth- 
erwise it is treated as a single pixel object and we assign a 
single temperature to the white dwarf. For the white dwarf 
we used pass band response functions of the UBVR filters. 

For reference of spatial structures in the disc, we use ra- 
dius and azimuth. The radius is used in units of the distance 
between the white dwarf and the inner Lagrangian point L\ , 
i.e. the point L\ has a radius of 1R L1 . The second coordinate 
is the azimuth, the angle as seen from the white dwarf and 
counting from —180° to 180°. Azimuth points towards the 
secondary, the leading lune of the disc has positive and the 
following lune negative azimuth angles. 



3.3 The distance 

Before reconstructions of physical parameters can be made 
reliably, the distance to the system must be known. Up to 
now, no reliable distance could be determined, because the 
secondary was not detected, not even in the infrared (Berri- 
man et al. 1986). The distance estimate of d — 90 to 150 pc 
exist with the assumption that the disc is optically thick 
(Watts et al. 1986). Allowing for optically thin emission in 
the disc, Berriman et al. find slightly larger distance esti- 
mates of between 120 and 170 pc. Line emission present in 
the quiescent state (Baptista et al. 1998) suggests, that the 
disc is at least partially optically thin (e.g. in the outer re- 
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Temperature Distribution 
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Figure 5. The reconstructed temperature 

gions, Watts et al.). A rough distance estimate of 184 pc, can 
be made using the white dwarf flux (Catalan 1998, private 
communication) . 

Vrielmann (2001) and Vrielmann, Hessman, Home 
(2002) describe the procedure how to arrive at a PPEM dis- 
tance estimate and successful applications of this method. 
Our method is based on the fact that at the best distance 
the PPEM reconstructions are smoothest, expressed by a 
maximal entropy. Fig. ^ shows the entropy as a function of 
the distance which peaks at 146 pc. For our further study 
we use the trial distance 145 pc which will give essentially 
identical results. 

The exact value of the final \ 2 of the fit to the observed 
data depends on the data set used, e.g. the amount of flick- 
ering still present after averaging or the applicability of the 
chosen spectral model. No standard value can be given for 
this final \ 2 ■ It is rather by choice of the analyst, judging 
the amount of structure in the reconstructions, the perfor- 
mance of the mapping algorithm and the quality of the fits 
and is therefore a rather subjective measure. It is therefore 
usual not to fit to a x 2 of 1, but to a somewhat higher value, 
allowing especially for deviations of the real spectral char- 
acteristic from the chosen spectral model. 

It is difficult to determine an exact error for the dis- 
tance, because of the uncertainty in difference between the 
spectral model used and the true spectral characteristic. In 
a future study we will attempt to quantify this error using 
more realistic spectra, such as calculated by Hubeny (1991). 
For the time being we estimate the PPEM distance error 
(with the assumption of a correct spectral model and good 
representation of the observation) generously to be about 



Surface Density Distribution 
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(left) and surface density (right) distributions. 

20 pc according to some limited testing with small changes 

in the constant part of the ephemeris (i.e. <pojQ correspond- 
ing to about twice the best Hipparcos distance estimates for 
objects with parallaxes less than 7 mas. 

The inclination angle i = 83.3°, the mass ratio q = 0.19, 
the white dwarf mass M w d = O.78.M0 and radius R w d = 
0.0244i?Q as well as the scale parameter R l1 /Rq — 0.42 
were taken from Baptista et al. (1998). 



3.4 The averaged light curves and final fits 

Fig. |^ shows the averaged quiescence light curves together 
with the final fit (\ 2 = 2.5) for a distance of 145 pc. While 
the overall profiles are reproduced, the fits deviate in detail 
from the observed light curves. This is mainly due to flares 
and flickering, e.g. at phases —0.135 and 0.13. The U-band 
and B-band light curves are fitted best, which means the 
Balmer jump is quite well reproduced. The residuals of the 
R-band light curve are worst, most likely due to the lack of 
line emission, particularly Ha, in our model spectrum. This 
agrees with the negative residuals at mid-eclipse and the too 
shallow model light curve. 

Further deviations, e.g. in disc ingress and egress, are 
probably due to the simplicity of the spectral model used. 



t An error in ipo shifts the white dwarf from the disc centre 
into the disc, introduces artefacts in the reconstructed map and 
thereby influences the entropy of the map as a function of trial 
distance. 
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Figure 6. Left: The azimuthally averaged intensity distributions I(T, S) (open circles) and Ibb(T) (filled triangles) for comparison with 
error bars representing the scatter in each annulus. The dashed line at 0.53 Rli indicates the size of the disc. Right: The ratio It,s/Ibb 
as a grey-scale plot. The solid line represents the Roche-lobe, the dashed circular line the size of the disc and the dotted lines ballistic 
stream trajectories for mass ratios 0.19 ±50%. The secondary Roche lobe is at the bottom, outside the plotted region. See text for further 
explanation. 



Keeping in mind that all light curves are fitted simulta- 
neously, we would be rather surprised that such a simple 
spectral model fits the data so well. A future application of 
a more sophisticated, more realistic model would probably 
give a better reproduction of the observed light curves. 



3.5 The reconstructions 

The reconstructed temperature and surface density distribu- 
tions are displayed in Fig. |^. In the centre the temperature 
reaches values up to nearly 15 000 K. It is fairly constant 
in the intermediate regions with a characteristic separation 
into a hot area with T ~ 10 500 K and a cooler region with 
T ~ 7500 K. Towards the disc edge it drops linearly. The 
surface density rises from the central parts to the disc edge 
from values around 0.02 to 0.03 g cm' 2 . 

As pointed out in VHH, it is essential to study the be- 
haviour of the spectral model in the given parameter space, 
before one can interpret the reconstructed maps. The easi- 
est disc parameter to estimate is the disc size. We just need 
to look at the intensity distribution derived from the T-E 
distributions, especially in the red filter. Since the R-band 
light curve is fitted worst, however, we would rather con- 
sider the V-band light curve in the present case. At a radius 
of 0.56-Rli the V-band intensity drops to about 1.5% of the 
maximum (i.e. the white dwarf intensity) and the distri- 
butions are to smoothed out due to the maximum entropy 
constraint. We therefore expect the disc edge to be between 
0.5077?^ and 0.567^1, i.e. at 0.53 ± 0.037^1. This disc ra- 
dius is compatible with being smaller than the 3:1 resonance 
radius R 3:1 = 0.68R L1 (for q = 0.19, Warner 1995), how- 



ever, significantly smaller than Steeghs et al.'s (2001) radius 
of 1.1R L1 (probably truncated due to tidal interactions to 
somewhere below 1R L1 ). Our disc radius describes the disc 
emitting continuum emission, while Steeghs et al.'s disc the 
line emission. It is possible that the chromosphere of the ac- 
cretion disc extends to larger radii than 0.56i?i,i, however, 
there is no indication (e.g. precession of the disc) that the 
quiescent disc is larger than 0.68R L i- We rather suggest that 
the disc as observed by Steeghs et al. was affected by the 
immediately preceeding outburst. 

When we compare the disc intensity I(T, E) (averaged 
over azimuth and wavelength) with the black body inten- 
sity derived from the temperature distibution alone Ibb{T) 
(Fig. ^) , we see in which parts the disc resemble black body 
emitters. In the inner part of the disc the two distributions 
are nearly identical and the ratio It,t./Ibb is very close to 
unity, i.e. the disc is nearly optically thick in the central re- 
gion up to a radius of about 0.15-Rli- Towards the edge the 
disc become more and more optically thin or the tempera- 
ture inversion in the chromosphere becomes more and more 
dominant. 

The grey-scale plot in Fig. ^ also shows that the nearly 
optically thick region extends towards the region where the 
bright spot is expected, i.e. where the accretion stream from 
the secondary hits the accretion disc. In this impact region 
the disc material will have a higher density compared with 
regions with the same radius but different azimuth. However, 
the bright region extends also towards smaller azimuths, i.e. 
spreads against the flow of matter in the disc. 

For a more quantitative analysis of the derived param- 
eters, we plot the averaged parameter distributions in Fig. [?] 
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Figure 7. Averaged temperature (left) and surface density (right) distributions. The error bars indicate the reliability of the reconstructed 
parameter values. At small radii the disc is opticlly thick and we can only give lower limits for the surface density, indicated by triangles. 



with error bars according to a study of the spectral model. 
The error bars are derived as follows: Similar to the study 
in VHH, we calculated the spectrum for the averaged pa- 
rameter pair (2\E) at a given radius, calculated spectra for 
a number of parameter pairs (T + AT, E + AE) around that 
parameter pair and determined the deviation of the spectra 
in terms of the x 2 . The error bars include the temperature 
and surface density range for a )( 2 of 3 for the given pa- 
rameter pair. Small error bars indicate a high sensitivity of 
the spectrum to the parameter, while large error bars mean 
only little variation in the spectrum in the given parameter 
range, e.g. in the optically thick range the surface density 
would have an infinitely large error bar. 

The error bars given in the averaged distributions 
(Fig. are good estimates for the error bars that should be 
attached to the values in the original distribution (Fig. 
They are omitted in Fig. [| mainly for clarity of the plot, 
but partly also because it would suggest an exact estimate 
of the error bars. 

We split the intermediate part of the disc (radii 0.2 to 
0.5R L i) into an "upper" and a "lower" branch. The upper 
branch corresponds to the bright spot region, where the disc 
is nearly optically thick, the lower to the other parts of the 
disc. 

The temperature is within the disc radius very well con- 
strained. In contrast, the surface density values are partly 
much less well defined. The bright spot region has large error 
bars, because it is nearly optically thick. At small radii the 



disc is optically thick and therefore we can only give lower 
limits for the surface density E. However, the true values of 
E can be very different from the given lower limits. 

The temperature drops with radius as expected. It 
reaches about 13 000 K in optically thick, inner part of the 
disc and drops to about 6 000 K at the disc edge. The opti- 
cally thick bright region has temperatures around 10 500 K 
while the other regions with identical radius have temper- 
atures between 9 000 K and 6 500 K. These temperatures 
are significantly lower than those quoted by Berriman et al. 
(1986) of 20 000 K. 

The surface density most likely increases with radius (as 
suggested by Fig. |H|) or could be roughly constant, consider- 
ing the partly large error bars in Fig. ^| (However, there is a 
clear distinction between the bright spot region and the rest 
of the disc.) Such a E(r) distribution is compatible with the- 
oretical calculations by Ludwig, Meyer-Hofmeister & Ritter 
(1994) and Cannizzo, Gosh & Wheeler (1982), but contra- 
dictory to Meyer & Meyer-Hofmeister's (1982) calculations. 
Since the critical surface density distributions within the 
disc instability model show an increase of E cri t (e.g. Can- 
nizzo & Wheeler 1984) we would rather expect that the ac- 
tual E(r) distribution also increases with radius. However, 
as we see later (Section 4.4), these surface densities probably 
describe only the chromosphere of the disc. 
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Table 4. The reconstructed fluxes (in mjy) for the uneclipscd 
component F un as well as a M4.5 and M5.5 main sequence star, 
F M4.5 and F M5.5' respectively. 



Filter 


Fun 


F M4.5 


F M5.5 


U 


0.59 






B 


0.44 






V 


0.00 


0.09 


0.01 


R 


0.66 


0.30 


0.07 



3.6 The white dwarf 

Using white dwarf spectra, the white dwarf temperature was 
reconstructed to 19 600 K using log g=8 white dwarf spectra. 
For the geometry we used only the "upper" half of the white 
dwarf, i.e. 0.25 + (90° - 83.3°)/360° = 0.269 of the total 
spherical surface. In the case that the disc is disrupted in the 
central part (i.e. if V2051 Oph is an intermediate polar) we 
would see more of the whole white dwarf and reconstruct a 
lower temperature. In case the white dwarf were fully visible 
the temperature would drop to 15 000 K. On the other hand, 
an accretion curtain may also hide part of the white dwarf 
surface, making it difficult to determine an exact value for 
the white dwarf temperature. 

We estimated the error of the white dwarf temperature 
Twd by checking how the light curve fit changes with differ- 
ent values of T m d if we keep the disc parameters as recon- 
structed. This gives us an error of about 1 000 K. However, 
the error in T w d is linked to the error in the distance. Taking 
this into account, the total error is rather 2 000 K. 

The geometry we used does not specifically set a loca- 
tion or size of a boundary layer. However, if such a bound- 
ary layer is present and emitting in the observed passbands, 
the emission should either be reconstructed in the innermost 
ring around the white dwarf or be hidden in the white dwarf 
emission. That we do not see a specific boundary layer emis- 
sion in the innermost ring of the disc either means that the 
white dwarf temperature is ov eres timated or that there is 
no boundary layer (see Section 4.6). 

Steeghs et al. (2001) and Catalan et al. (1998) both 
find a white dwarf temperature of 15 000 K. However, since 
both do not state which distance value they used a direct 
comparison to our value is not possible. Assuming they used 
a similar distance, this would mean that we might see part of 
the "lower hemisphere". However, this is doubtfull, since the 
disc is optically thick near the white dwarf and should occult 
the "lower" half of the compact source. If the boundary layer 
were hidden in the white dwarf emission, Steeghs et al. and 
Catalan et al. should have faced the same problem. 

The reconstructed value is close to typical values of 
white dwarf temperatures, but somewhat higher than in 
the similar objects HT Cas, OY Car and Z Cha of around 
15 000 K (Gansicke & Koester 1999 and references therein). 
This might be connected to the somewhat different system 
parameters as shown in Tab. ^ 



any spectral model. It is simply any uneclipsed flux in the 
given filter and free for any interpretation. 

Since the uneclipsed spectrum might be generaly under- 
estimated, the B-flux overestimated (Vrielmann & Baptista 
2002) and since the errors in the uneclipsed component par- 
ticularly in the V and R-band filters are rather large (Vriel- 
mann et al. 2002), we cannot derive much physics out of this 
component. All we can say is that the uneclipsed component 
is compatible with a M5 main sequence star (using Baptista 
et al.'s mass and radius and Kirkpatrick & McCarthy 1994 
tables for main sequence stars, see Table ^) plus an opti- 
cally thin component emitting particularly in the U-band 
and R-band, e.g. a chromosphere and/or disc wind. 

Evidence for a disc wind comes from the fact that the 
blueshifted Doppler emission component of all emission lines 
is always brighter than the redshifted component (Kaitchuck 
et al. 1994, Honeycutt, Kaitchuck & Schlegel 1987). Even 
though Steeghs et al. (2001) claim to find a reversal of this 
asymmetry for Ha (i.e. the red peak is stronger) we can- 
not see this in their gray-scale plot: their Fig. 1 shows the 
blue wings stronger at all plotted phases (-0.3 to 0.1) except 
during part of the eclipse. Furthermore, the eclipse is shorter 
and shallower in the blue peak than in the red peak, indicat- 
ing that the origin of the blue emission is spatially focused 
and only partly eclipsed, e.g. like a collimated wind (or jet) 
above the disc. And although in Warner & O'Donoghue's 
(1987) spectrum the red peak is higher, the blue peak is 
slightly broader, indicating that the additional emission in 
the blue wings has rather high velocities. 

3.8 Derived parameters 

From the reconstructed parameter distributions we can cal- 
culate several further parameters, like the effective temper- 
ature, the disc scale height and the viscosity. These param- 
eters are more descriptive of the disc than the gas temper- 
ature and surface density. The disc size and effective tem- 
perature can in particular be determined more accurately 
than the originally reconstructed ones, because they rather 
depend on the intensity distributions corresponding to the 
original reconstructions than on the reconstructed parame- 
ters themselves. However, in general the derived values are 
only as certain as the originally reconstructed values on 
which they are based and one has to be careful in interpret- 
ing the resulting derived parameter values. One must also 
keep in mind, that these parameters might only des crib e 
the hot chromosphere, not the entire disc (see Section 4.4). 
Furthermore, it is difficult estimating an error of the derived 
parameters without extremely time consuming Monte-Carlo 
type calculations. 



3.8.1 The effective temperature and mass accreion rate 

From the temperature and surface density distributions we 
can calculate the spatially resolved effective temperature T e a 
distribution in the disc (Fig. ^). T c ff is calculated from the 
intensity 7„(T, £) as 



3.7 The uneclipsed component 

Apart from the disc parameters we also reconstructed an 
uneclipsed fraction of the light curve (Table ^). For the re- 
construction of this component, we do not need to assume 



ty / I V {T, E) dv — <jT cS 



(3) 



with a the Stefan-Boltzmann constant. In practice, I v is 
summed numerically over a large frequency range v. Equa- 
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Effective Temperature Distribution 




Radius (R L1 ) distance (R L1 ) 

Figure 8. Left: The radial effective temperature distribution in quiescence. The underlying dotted lines are theoretical steady state 
temperature distributions for log Al = 13 to 18, the one for M = 10 16 gs — 1 is drawn solid for reference. The dashed-dotted lines indicates 
the critical effective temperature T c g ^ according to Ludwig, Meyer-Hofmeister & Ritter's (1994) critical mass accretion rate M\. Right: 
The spatial distribution of the effective temperature as a grey-scale plot. The solid line represents the Roche-lobe of the primary, the 
dashed line (in both plots) the disc size and the dotted lines ballistic accretion stream lines for mass ratios 0.19 ± 50%. The secondary 
Roche lobe is at the bottom, outside the plotted region. 



tion shows that T e s does not directly depend on the recon- 
structed values T and E, but rather on the intensity derived 
from them. Since the intensity distribution directly trans- 
lates into the observed fluxes which are as good as the fit to 
the light curves suggests, the effective temperature is more 
reliable than the surface density values and rather compara- 
ble to the certainty of the reconstructed temperature values. 

The main part of the disc has a mass accretion rate 
of between M = 1 x 10 15 gs" 1 to M = 1 x 10 17 gs" 1 
(1.6 x lCT^/Vf© yr -1 to 1.6 x 10~ 9 M Q yr" 1 ). However, the 
distribution is split into two parts of higher and lower tem- 
peratures. As the grey-scale plot in Fig. ^| shows, the higher 
temperatures are associated with the bright spot region, or 
more general to the region between the two stars, i.e. the 
same region that is nearly optically thick. The remainder of 
the disc is cooler. 

This situation is understandable in terms of the mass 
accretion rate involved. At the bright spot the mass transfer 
is highest, reaching 1 x 10 17 gs -1 , because of the impact of 
the gas stream from the secondary. In the remaining disc, the 
mass accretion rate is determined solely by the viscosity. In 
the region 180° away from the gas stream the mass accretion 
rate drops to values around 3 x 10 gs - . 

Fig. |^ also shows the theoretical critical effective tem- 
perature T e ff t A(R) derived from the critical mass accretion 



rate Ma(R) as given by Ludwig et al. (1994), indic atin g 
whether the disc should undergo outbursts (see Section 4.5). 



3.8.2 The disc mass 

The disc mass can be estimated from the surface density, 
however, only with a certainty that the reconstructed surface 
density values allow. For the disc within a radius of 0.53-Rli 
we derive a mass of 7.7 x 10 19 g or 3.9 x 1O" 14 A1 . This 
value has to be treated with some caution. The disc is nearly 
optically thick in the centre, i.e. for that region we rather 
derive a minimum value of E and therefore of the disc mass. 
In the outer parts of the disc the surface density might vary 
by up to a third of the reconstructed value for E. The error 
in the disc mass therefore must be set to roughly 50% of the 
derived value. 

According to Frank, King & Raine (1992) the mass in 
the disc at any one time is M disc < 10~ 10 M Q a~ 4/5 A^ 1 °, 

where Mie is the mass accretion rate in units of 10 16 gs -1 
and a the viscosity coefficient determined by u = ac 3 H, 
where v is the viscosity, c 3 the sound speed and H disc 
thickness. So, unless a reaches unreasonably high values of 
the order of 10 4 or the underlying disc is extremely massive, 
Frank et al.'s condition is well fulfilled. 
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Disc Height Distribution 




Figure 9. The radial distribution of the disc scale height. The 
circle at a disc radius of ORli represents the radius of the white 
dwarf and illustrates how geometrically thin the disc is. The scale 
of the radius (x-axis) would have to be multiplied with a factor 



of 2.9 X 10 10 cm Rl 
scale height. 



to be translated into the same units as the 



3.8.3 The disc scale height 

We can estimate an equilibrium scale height of the disc (or 
rather its chromosphere, Fig. ^). The value of H/R ranges 
around 0.01, i.e. the opening angle of the disc is about 0?6 
with a scatter of 0? 1 and the disc has a scale height of only 
1/4 white dwarf radii at the disc radius. The error on the 
values of H are fairly small (less than 1%), as this parameter 
only depends on the square-root of the temperature. This 
validates our assumption of a geometrically thin disc. Up to 
a radius of 0.4i? L i the disc has a slight concave shape with 
the opening angle ranging between 0?4 and 0°85. 

However, this disc scale height is calculated with the 
assumption of an iso-thermal vertical (z-)structure of the 
disc. Realistic calculations of the disc photosphere lead (for 
mass accretion rates estimated in this paper) to disc opening 
angles of about 3 degrees in addition to the underlying disc 
with unknown, but most likely non- negligible opening angle. 
In an as high-inclination system as V2051 Oph this should 
lead to significant front-back asymmetries in the disc emis- 
sion (however, opposite to what we observe, see Section 4.3). 



3.8.4 The disc viscosity 
Using the formula 



-Q K aPH = / F v dv 



(4) 



where Qk = \J GMw / R is the Keplerian angular velocity, 
P the gas pressure (calculated according to the perfect gas 
law: T = fimnP/ pk, with T the temperature, \i the ioniza- 
tion coefficient, ran the mass of the hydrogen atom, p the 
density and k the Boltzmann's constant), H again the disc 
scale height and F u the flux at frequency v, we can deter- 



mine the viscosity parameter a (Shakura & Sunyaev 1973) 
throughout the disc. Typical values of a in a quiescent disc 
should lie around 0.01 (e.g. Smak 1984). 

The derived values are very high (between 50 and 900) 
and incompatible with current accretion disc models. A very 
similar situation arises in the dwarf nova HT Cas (Wood, 
Home & Vennes 1992, Vrielmann 1997, Vrielmann et al. 
2002). We see a need to explain these high values. 

We suggest that the structure of the disc is such, that we 
do not see the entire disc. For example, in the case that part 
of the disc is obscured or outshined by a hot chromosphere, 
the gas density is higher than calculated and a would be 
lower. In Section 1.4 we propose a model of the disc, which 
could - in principle - explain the observations. However, 
only concrete theoretical calculations of such a scenario can 
show if such a model is indeed applicable. 



4 SUMMARY & DISCUSSION 
4.1 The distance 

We derived a new distance to V2051 Oph of 146T20 pc using 
the Physical Parameter Eclipse Mapping (PPEM) method, 
assuming that the true distance would lead to the smoothest 
reconstruction. Our result is compatible with those of Watts 
et al. (1986) and Berriman et al.'s (1986). This agreement 
and a few others others (Vrielmann 2001) support the ap- 
plicability of PPEM to derive distances even when we use a 
very simple spectral model for the disc emission. However, 
the method certainly has to be applied to a large number 
of systems and compared to independent distance estimates 
before it can be considered established. 



4.2 The radial light distribution of the accretion 
disc 

We applied the PPEM to a quiescent UBVR data set. The 
disc appears nearly optically thick inside a radius of about 
0.2R L1 and in the region where the bright spot is expected 
(up to a radius of about 0AR L1 ), while outside these regions 
the disc is optically thin or shows an increasing dominance 
of a temperature inversion in the disc chromosphere. This 
agrees with Tylenda (1981) and Watts et al. (1986) in gen- 
eral, although not in detail. 

Tylenda's (1981) calculations show that a disc with a 
mass accretion rate of 10 16 gs _1 is optically thin in the outer 
parts at radii larger than 1.2 x 10 10 cm (for V2051 Oph this 
corresponds to 0.42i?i,i). Since the accretion disc of V2051 
Oph has a radius of about 1.55 x 10 10 cm (~ 0.53J?i,i) this 
transition zone lies well within the disc radius, however, is 
incompatible with ours. 

Using the equivalent width (EW) of the Balmer emis- 
sion lines during an orbit, Watts et al. found the radius of the 
transition between the dominance of optically thin over thick 
regions to be at a radius of 0.16±0.02a (= 0.24±0.03fln). 
Their value lies inbetween our radius of the optically thick 
part of the disc (0.2 and 0.4J? L i), because their method im- 
plies axisymmetry of the accretion disc emission. 

The fact that emission lines of V2051 Oph have broad 
wings that originate at radii R < 0.2R L1 (if the disc is Kep- 
lerian, see also also Section 4^6) must mean that they orig- 
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inat e fro m a (chromospheric) layer above the disc (see Sec- 
tion 4.4). 



4.3 The bright region in the accretion disc 

Why does the disc show the bright region between —30° and 
+60° (Fig. ^) as well as a high temperature ridge from the 
white dwarf to the edge of the disc at an azimuth of about 
-10°? 

This bright region is associated with enhanced emission 
and therefore not explicable as a back-front assymetry as 
expected for a disc with a significant, underestimated (or 
neglected) opening angle. The latter would lead to an en- 
hancement in the "back" part of the disc, i.e. at azimuths 
around ±180°. 

The bright spot region appears to be nearly optically 
thick in quiescence and slightly preceding the theoretical 
position in the direction of motion in the disc. However, one 
would rather expect the bright spot region to extend from 
the impact area into the direction of motion in the disc. Even 
though this bright region might be only temporary and is 
only derived using observation near the eclipse, we discuss 
its possible origin. 

We carefully checked if the reconstructed location of the 
bright spot region is caused by an errornous ephemeris, but 
must conclude that this cannot be the case. This means, if 
the feature is real, there must be a physical reason for this 
unexpected disc structure. Since the secondary plays only a 
minor role in V2051 Oph, as its contribution is small com- 
pared to the disc, it cannot be responsible for illumination 
of the disc. 

The high temperature ridge may indicate that the disc 
is warped and the ridge at azimuth —10° maximally illumi- 
nated by the white dwarf. The azimuths +170° would then 
not be illuminated by the white dwarf. Since the data we 
used for the PPEM analysis cover only the phases ±0.15 we 
never see the outer regions of the disc around azimuth angles 
±180° leading to low temperatures in our reconstructions. 
This warp would, however, have to be fixed in the rotating 
frame to be so prominent in the maps derived from averaged 
light curves. But calculations by Schandl & Meyer (1994) 
who propose coronal winds as the cause of the warp, show 
a precession of the disc, making the warp model unlikely. 

Another possible scenario is a bulge in the disc similar as 
the one proposed by Buckley & Tuohy (1989) at the location 
of the impact region. This bulge could illuminate the disc 
area around azimuth angle 0. How the high temperature 
ridge can be formed is unclear and it would be interesting 
to analyse multi-colour observations of other epochs to see 
if this feature is stable. 

Kaitchuck et al. (1994) calculated Doppler images of 
the emission lines H/3 and He I A4471. They find that the 
leading side of the disc is brighter than the trailing side. 
The location of their bright area does not coincide with our 
bright region. However, they investigate the location of the 
line emission, which should be minimal in an optically thick 
area and could therefore correspond to the bulge itself or the 
region ahead of the bulge (i.e. at even larger azimuths) as 
proposed by Hoard et al. (1998) for the nova- like UU Aqr. 
Their disc model involves an emitting impact spot with a 
trail of optically thin material along the edge of the disc. 

The size and location of the bright region will be cor- 



related with the mass accretion rate M, since the latter 
determines the disc size. It would therefore be interesting 
to study V2051 Oph at various times during high and low 
states. For example, the variability in M may also be the 
reason, why Watts et al. (1986) could not detect the s-wave 
(produced by a bright spot), while it was clearly present in 
Honeycutt, Kaitchuck & Schlegel's (1987) data and might 
indicate a variability in the bulge. 

4.4 The vertical disc structure 

We propose that the disc itself must have a structure as fol- 
lows, and similar to the disc structure suggested for HT Cas 
(Vrielmann 2001, Vrielmann et al. 2002): The disc consists 
of a cool, optically thick layer with a hot chromosphere on 
both sides of the disc. This chromosphere is nearly optically 
thick at small radii, but optically thin for most of the disc. 
What we see of the disc is mainly the hot chromosphere and 
little from the cool layer underneath, therefore the temper- 
ature and in particular the surface density are descriptive 
only of the chromosphere. Such a model was first proposed 
by Wood et al. (1992), Menou (2002) has done some more 
recent theoretical work on it. 

This disc model resolves in the confusion about the op- 
tical depth of the disc and subsequently our derived values of 
a. On the one hand, the disc instability model (e.g. Ludwig 
et al. 1994, Faulkner et al. 1983, Smak 1982) predicts that 
the surface density values should lie between 10 and 100 g 
cm~ 2 for a's between 0.1 and 1 and even at larger E's for 
smaller a's. This means that the disc must have a massive, 
optically thick component. 

On the other hand we see line emission from the disc 
during quiescence, rise and decline of an outburst (e.g. Hes- 
sman et al. 1984) which indicates that there must be an 
optically thin component during these phases in the out- 
burst cycle. Our derived surface density values in the disc 
re-confirms that the optically thin source (or more correctly, 
the non-black body source) is associated with the disc. 

The sandwich model elegantly combines both the ob- 
servational findings (optically thin material) and the theo- 
retical predictions (optically thick disc) . It also explains why 
we find such high values of a: Since we see mainly the con- 
tribution from the hot chromosphere, the surface densities 
describe only the upper layers of the accretion disc. The true 
surface densities (E) must be much larger, of the order 10 to 
100 g cm' 2 . A simple estimate confirms that values around 
100 g cm -2 for £ lead to a's in the expected range of 0.01 
to 1. Thus, we have found a simple explanation, why we de- 
rive such high values for a without discarding the accretion 
disc models or our PPEM method: Using the surface den- 
sity of the chromosphere alone and assuming it to describe 
the whole accretion disc leads to high values of the viscosity 
parameter a. 

The proposed chromosphere must extend to such 
heights above the orbital plane that part of it is never 
eclipsed, i.e. near the white dwarf it must reach farther 
up t han about 0.2R L1 or O.O85J?0 or about 87?^,^ (see Sec- 
tion 3.7). This chromosphere should also be the site of the 



observed emission lines which would be produced by non- 
radiative energy dissipation, analogous to a stellar chromo- 
sphere (Home & Saar 1991). 

A more extensive explanation of the vertical disc struc- 
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ture is given in Vrielmann et al. (2001) for the similar system 
HT Cas. 

4.5 The mass accretion rate 

Similar to results of Watts et al. (1986), the quiescent mass 
accretion rate appears in large parts of the disc to be some- 
what higher than or around the critical values given by Lud- 
wig et al. 1994 or Watts et al. of ~ 4.6 x 10~ n A^q yr" 1 (for 
K\ = 111km s~l). The latter derived a similar accretion 
rate for V2051 Oph as we of ~ 1.3 x 10 _10 X Q yr -1 com- 
pared to our 1.6 x 10 _11 X Q yr^ 1 to 1.6 x 10" 9 Xq yr _1 . 
Fig. ^ shows that specifically the bright spot and inner, hot 
regions (R < 0.1R L1 ) have higher temperatures than the 
critical values T a s t A(R)- 

Since V2051 Oph shows outbursts, and probably very 
infrequently, the effective temperature should not lie so close 
to or even above the critical T c h : a{R), but rather far below 
(e.g. Warner 1995). Berriman et al. (1986) also found that 
V2051 Oph's disc is too hot to undergo the disc instability 
cycle. The only known outbursts around the time of our 
observations (June 1983) were the ones in April 1982 and 
in May 1984. V2051 Oph has infrequent observed outbursts 
and a typical outburst cycle of about a few years (VSNET 
observations). So the disc is unlikely to have been affected 
by a rise to or decline from outburst. A similar situation 
occurs in HT Cas (Vrielmann et al. 2002), a dwarf nova with 
reportedly infrequent outbursts. At the moment, we have 
no solution for this problem. Until another explanation is 
found, we must assume that either our mass accretion rates 
are errorneous or that it is sufficient if part of the disc has 
effective temperatures below the critical value. 

We would recommend a study of this object throughout 
a quiescent cycle in order to check variations in the mass ac- 
cretion rate, especially before the onset of an outburst and 
during low states (B ~ 16.2, Baptista et al. (1998), com- 
pared to B ~ 15 (WO). Interestingly, HT Cas, a similar 
dwarf nova, also shows variations during quiescence (Robert- 
son & Honeycutt 1996) with fluctuations of about 1.5 mag. 

Furthermore, it would be very desirable to get an in- 
dependent distance estimate. If the distance were signif- 
icantly smaller than our estimate, the effective tempera- 
tures and mass accretion rates would be reconstructed to 
smaller values. However, this will not solve the problem, be- 
cause an (unlikely) reduction of a factor two in distance will 
only lead to a small reduction in effective temperatures to 
0.5 1 ^ 4 T c ff = 0.84T e ff and therefore to a small change in M. 

4.6 Is V2051 Oph an intermediate polar? 

While it is clear that V2051 Oph shows typical features of a 
SU UMa system (outbursts and super outbursts), it is not 
clear, whether the white dwarf has a small magnetic field. 

An intermediate polar should show a disrupted inner 
disc in the PPEM maps, unless the inner disc radius is too 
small. We do not see any indication for a hole, except for 
a flattening of the effective temperature towards the white 
dwarf. Only a very tiny hole could be hidden due to the 
MEM smearing, as tests of the PPEM method show. How- 
ever, estimating the line width at phase 0.51 of Ha in WO's 
Fig. 10, we yield ~ 77A, i.e. the disc material at the in- 
ner disc edge has a velocity of about 1777 kms -1 . For a 



Keplerian disc this corresponds to an inner disc radius of 
3.3 x 10 9 cm = 0.11-Rli. More recent, yet unpublished data 
of V2051 Oph suggest line wings extending to about 2530 
kms -1 , corresponding to an inner disc radius of 2.3 white 
dwarf radii. Quiescence data of Watts et al. (1986) also show 
line wings of H/3 and H7 extending to maximally 2500 kms -1 
(corresponding to R in — 0.06R L1 = 2AR m d). 

Steeghs et al. (2001) find much broader line wings. Both 
spectra sets were taken during decline from outburst which 
means that the emission lines should have been affected in 
the same way. Since the intensity in the outer line wings 
is extremely low (< 3% for v < —2500 kms -1 ; red wing 
affected by He line, Fig. 3 of Steeghs et al.), broad shal- 
low wings could - in principle - be lost in the noise of the 
above mentioned spectra, however, no such line profiles as 
observed by Steeghs et al. can be hidden in those spectra. 
Only the blue wing of H/3 can be measured reliably, which 
will be enhanced by the disc wind, a possible explanation 
for Steeghs et al.'s broad line profile. 

The fact that we do not see a hole in the reconstructed 
disc means either, there is none, it is extremely small, or it 
could be hidden by accretion curtains. 

X-ray observations (Holcomb, Caillault & Patterson 
1994, van Teeseling, Beuermann & Verbunt 1996) and lack 
of circular polarisation (Cropper 1986) data do not contra- 
dict the intermediate polar model with a very weak magnetic 
field of 1 MG or less, as expected for intermediate polars. 

Thus, if V2051 Oph is indeed an intermediate polar, its 
white dwarf must have a very weak magnetic field which 
allows to disc to be disrupted with an inner disc radius of 
Rin < 0.1R L1 = 4:.2R w d- At present we cannot distinguish 
between such a low field intermediate polar and a dwarf 
nova. However, the variable hump and irregular occurance 
of strong flares as observed by WO and Berriman et al. 1986 
might indicate a flaring accretion curtain which lights up 
during inhomogeneous accretion. Warner & Woudt (2002) 
propose a model for V2051 Oph in which the spinup of the 
magnetic field lines of a small field on the white dwarf tem- 
porarily creates such a small hole of a a few tenths larger 
than the white dwarf radius. Their model applies to a disc 
on decline from outburst. During quiescence they expect the 
hole to be even slightly larger. If this hole is present and per- 
sisting long enough, accretion curtains can develop, leading 
to the observed flares. 

Simultaneous multi-colour observations at different 
epochs and high resolution spectroscopy are very desirable 
in order to solve the mystery of V2051 Oph's nature. 

4.7 Comparison to similar objects 

V2051 Oph has a similar orbital period as OY Car, Z Cha 
and HT Cas and since all of them are SU UMa system, show- 
ing super outbursts inbetween their normal ones, one would 
expect to find some similarities among these systems. How- 
ever, the light curve of V2051 Oph seems to vary a lot more 
than that of the other three systems, especially in terms of 
the presence and location of the (orbital) hump, flickering 
and flaring. 

The eclipse light curve of OY Car (Vogt et al. 1981, 
Wood et al. 1989, Rutten et al. 1992, Home et al. 1994) is 
remarkably stable over more than a decade and even close 
to a super outburst (Hessman et al. 1992). A similar consis- 
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tency in the appearance is shown by HT Cas (Berriman et 
al. 1987, Home et al. 1991, Robertson & Honeycutt 1996) 
and Z Cha (Wood et al. 1986, Marsh et al. 1987, Robinson 
et al. 1995). OY Car and Z Cha even show text book light 
curves of eclipsing dwarf novae in which the orbital hump 
and the ingress and egress phases of the white dwarf and the 
bright spot can be very well distinguished and explained in 
the standard Roche model. HT Cas, on the contrary, has 
never been observed to show an orbital hump. 

For comparison, we listed a few system parameters of 
these four objects in Tab. []. While at first glance all pa- 
rameters seem to be similar, closer inspection reveals that 
the secondary in V2051 Oph is about twice as massive as 
the secondaries in the other systems. Furthermore, also the 
mass accretion rate M into the disc of V2051 Oph is larger 
by a factor of twenty or more than in the other systems 
(in OY Car and Z Cha M was determined by Eclipse Map- 
ping (Wood et al. 1989 using d=100 pc, Home et al. 1994; 
Wood et al. 1986) and in HT Cas with the PPEM method 
(Vrielmann et al. 2002)). 

The larger amount of variation of V2051 Oph could be 
induced by the larger mass accretion rate, being close to 
or above the theoretical limit. Even if the disc were rather 
patchy with a covering factor of maybe 50%, as proposed 
for HT Cas (Vrielmann et al. 2002), the mass accretion rate 
were still a factor of about ten or more larger than that of 
the other three dwarf novae. The larger mass of the sec- 
ondary in V2051 Oph, could be the reason for the higher 
mass accretion rate, however, it is not fully understood, how 
Ai is related to other system parameters. Furthermore, it is 
not clear, how the larger M alone can cause the occational 
occurence of a hump. We therefore, cannot exclude the low- 
field intermediate polar model. 

4.8 Concluding remark 

Since V2051 Oph is usually too faint for amateur as- 
tronomers, monitoring of this object is difficult, but partic- 
ularly interesting. Especially also the evolution of the disc 
during normal and low quiescent states and the percentage 
of time V2051 Oph spends in either state would be of im- 
mense interest. 
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Physical Parameter Eclipse Mapping of the quiescent disc in V2051 Ophiuchi 



Table 5. System parameters of V2051 Oph and similar SU UMa type dwarf novae in comparison. 



Parameter 


V2051 Oph 


HT Cas 


OY Car 


Z Cha 


P [h] 


1.50 


1.77 


1.51 


1.79 


i 


83?3 ± 1?4 


81?0 ± 1?0 


83?3 ± 0?2 


81?7 ± 0?13 


q 


0.19 ± 0.03 


0.15 ± 0.03 


0.102 ± 0.003 


0.1495 ± 0.0035 


M wd [M Q ] 


0.78 ± 0.06 


0.61 ± 0.04 


0.685 ± 0.011 


0.54 ± 0.01 


Mr [M e ] 


0.15 ± 0.03 


0.09 ± 0.02 


0.070 ± 0.002 


0.081 ± 0.003 


M [gs" 1 ] 


1 x 10 17 


5 x 10 15 


4 x 10 15 


4 x 10 15 



M: Mass accretion rate at the disc radius (disc edge) or into the bright spot 

References: V2051 Oph: Baptista et al. (1998), HT Cas: Home et al. (1991), Vriclmann et al. (2002, M), 
OY Car: Wood et al. (1989), Home et al. (1994), Z Cha: Wood et al. (1986) 



Vriclmann, S., Baptista, R., 2002, submitted to AN 
Vrielmann, S., Home, K., Hcssman, F.V. 1999, MNRAS, 306, 766 
(VHH) 

Vrielmann, S., Hessman, F.V., Home, K. 2002, MNRAS, in press 
Vogt, N., Schoembs, R., Krzcminski, W., Pedersen, H. 1981, 

A&A, 94, L29 
Tylenda, R. 1981, Acta Astron., 31, 127 

Warner, B. 1995, "Cataclysmic variable stars", Cambridge Astro- 
physics Series, Cambridge University Press, p. 207 
Warner, B. 1996, Ap&SS, 241, 263 
Warner, B., Woudt, P.A. 2002, submitted to MNRAS 
Warner, B., Cropper, M. 1983, MNRAS, 203, 909 
Warner, B., O'Donoghue, D. 1987, MNRAS, 224, 733 (WO) 
Watts, D.J., Bailey, J., Hill, P.W., Greenhill, J.G., McCowage, 

C.,Carty, T. 1986, A&A, 154, 197 
Wood J.H., Home, K., Berriman, C, Wade, R., O'Donoghue, D., 

Warner, B. 1986, MNRAS, 219, 629 
Wood J.H., Home, K., Berriman, C, Wade, R. 1989, ApJ, 341, 
974 

Wood J.H., Home K., Vennes S., 1992, ApJ, 385, 294 



© 0000 RAS, MNRAS 000, 000-000 



